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In this work, we present new soliton solutions for thick branes in 4+ 1 dimensions. In particular, 
we consider brane models based on the sine-Gordon (SG), and scalar fields, wliich have 
broken Z 2 symmetry in some cases, and are responsible for supporting and stabilizing the thick 
branes. The origin of the symmetry breaking in these models resides in the fact that the modified 
scalar field potential may have non-degenerate vacuua. These vacuua determine the cosmological 
constant on both sides of the brane. We also study the geodesic equations along the fifth dimension, 
in order to explore the particle motion in the neighbourhood of the brane. Furthermore, we examine 
the stability of the thick branes, by determining the sign of the term in the expansion of the 
potential for the resxilting Schrodinger-like equation, where w is the 5-dimensional coordinate. It 
turns out that the brane is stable, while there are unstable modes for certain ranges of the 
model parameters in the SG and branes. 
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I. INTRODUCTION 

Since there is no known fundamental principle requiring spacetime to be (3 + 1)— dimensional p][2], it has been 
suggested that our observable universe might be a (3+ 1)— dimensional brane in a higher dimensional space [3][3. In 
most models, there are one or more fiat 3-branes embedded discontinuously in the ambient geometry |5]. Moreover, 
ideas with two S-braues provide a very elegant description of the large hierarchy between the scales of weak and 
gravitational forces ^ [7| and contain massless modes wliich reproduce Newtonian gravity at large distances on the 
brane In recent years, particle physics extra-dimensional theories beyond the standard model have become a 
standard part of the array of ])henomenological models [51 [§]. Although there are still no experimental evidence 
supporting extra dimensions, due to various theoretical motivations, extra dimensional models continue to be widely 
considered in the literature [T^. In this context, it would be quite useful to have a set of simple and sufficiently 
general rules which would allow one to test new models 0. Most extra-dimensional models require the existence of 
scalar fields, for instance, to generate a domain-wall which localizes matter fields [3]. The scalar fields also serve to 
stabilize the size of the compact extra dimensions [SlSi niochfy the Randall-S undr uni warped-space 

(^02] to a smoothed-out version 0E]. or to cut off the extra dimeiLsion at a singularity [Jll^. In order to replace 
an infinitely thin brane with a thick one, a scalar field with soliton behavior is frequently invoked. The nonlinearity in 
the scalar field and in particular the existence of discrete vacuua in the self-interact ion of the scalar field lead to the 
appearance of a stable localized solution, which is a good motivation for building thick brane models {TO]* A general 
method for determining the lowest energy configuration lias been worked out in [m 1T3] . 

Recently, brane world models have also been considered in liigher order curvature gravity and in modified teleparallel 
gravity. For iiLStance, five-dimensional modified teleparallel gravity was considered in a brane scenario, where analytic 
domain walls were found to have a double-kink solution in the aftermath of the torsion of spacetime [15]. Furthermore, 
this model was extended by using a first-order formalism to find analytical solutions for models that include a scalar 
field with standard and generalized dynamics. In addition to this, it was found that the brane splits, as a result of 
the deviation from the standaixl model by controlling specific parameters [IS]. Modified gravity in five dimensional 
spacetime has also been analysed in the Palatini formalism. For instance, a thick Palatini /(i?) brane described by 
an anti-de Sitter warped geometry with a single extra dimension of infinite extent, sourced by a real scalar field was 
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stiidieil ill a perturbative scenario (T7]. Besides, tlie model of a domain wall {thick brane) in a non-compact AdS 
space time with only one extra dimeiLsion was further analysed in [T§1[^. The cla^isical tasts of General Relativity in 
thick branea were also studied by studying the motion of test particles in a thick brane scenario and the impact of the 
brane tliickness on the four-dimeiLsioiial path of massless particles was explored in |2n]. More speciBcaUy, by applying 
a confinement mechanism of massive tests particles in the domain wall, for instance, that simulates classically the 
trapping of the Dirac field in a domain wall, the influence of the brane thickness on the four-dimensional (4D) path of 
massless particles was analysed. A generalized version of the Randall-Siiiidrum II model with different cosmological 
constants on each sirle of a brane were also discussed, where specific configuratioiLS of a scalar field and its stability as a 
replacing factor of the singular brane were coiLsidered [^. Models of thick branes in noncompact five-dimensional bulk 
with different anti-de Sitter geometries on eacli side of the brane were explored [221[23]i cosmological applications 
of soliton-like thick branes have also been studied |23]. On the other hand, the existence of brane solutions were also 
considered as a result of a real scalar' field in the presence of five dimensional f{B) gravity (2S]. In addition to this, 
asymmetric thick branworld scenarios were studied, by changing the superpotential of the scalar field [2S1« 

It is widely practiced that in brane world scenarios Z 2 symmetry is assumed [ElEZlESli wliich is originally motivated 
from the Z 2 syninietry considered in M-theory [2Z|. Under this symmetry the bulk metric on the two sides of the 
brane should be the same [SHI* Moreover, under sucli a symmetry the empty bulk on either sides of the brane 
have the same negative cosmological constant and as a result they are AdS [2H]« Note that these conditions are 
satisfied in the Randall-Sandrum model. There are, however, brane models in whicli there is no Z 2 symmetry and the 
bulk is different on both sides of the brane (23. In the latter, the Friedmann equation for a pasitive brane tension 
sittiated between two bulk spacetimes that posses the same 5D cosmological constant, but which does not possess a 
Z 2 symmetry of the metric itself was derived, and the passible effects of dropping the Z 2 synmietry on the expansion 
of our Universe were examined. In some of these models, the cosmological constant differ on both sides of the brane 
m< where the effects of including a Gauss-Bonnet combination of higher-order curvature invariants in the bulk action 
are taken into accoimt. In fact, by considering brane world scenarios including the Gauss-Bonnet term, it was found 
that the cosmological dynamics have the same form as those in Rand all-S undr uni scenarios but with time-varying 
four-dimensional gravitational and cosmological constants |23. Motivated by such a possibility, we consider in this 
paper several models, namely, the sine-Gordon (SG), and brane models which have broken Z 2 symmetry in 
some c?Lses. What is meant in the paper ?ls Z 2 syimuetry is the symmetry with respect to the position of the brane 
(not the actual symmetry in the lagrangiau). The brane position is shifted from 2 = 0 for some models. As a result, 
we find that for the 0® and SG systems this symmetry is broken and the vacua on the two sides of the brane are not 
degenerate. In several cases, the Z 2 symmetry can be restored by a proper choice of model parameters. The origin of 
symmetry breaking in our models resides in the fact that the modified scalar field potential may have non-degenerate 
vacuua. These vacuua determine the cosmological constant on both sides of the brane. 

Relative to the stability issue, topological soli tons are known for their non-singulai' structure and a natmal local¬ 
ization mechanism which are highly stable. Zeldovich et ai [TTl E3 suggested that the soliton of the model is a 
rea.sonable source for the formation of domain walls. Vilenkin extended this idea to incorporate the general theory 
of relativity (TTISIIES* Since there is a clase similarity between domain walls in 3+1 dimensions and branes in 4+1 
dimensions, it is natural to think that the soliton idea might have something to do with the existence and stability of 
branes. Motivated by this idea, we consider soliton models for thick branes and extend some of the existing works. In 
particulai', in [25] the authors have studied the stability of the kink brane model in five dimensioiLS, as well as the 
and the inverted <p‘* potential. Furthermore, the properties of fermions coupled to the sine-Gordon brane model 
were inve.stigated. In [23. double kink-like solutions were considered and the stability of scalar, vector and tensor 
perturbations were discussed. 

Tliis paper is outlined in the following manner: In Section [T^ we briefly review the thick brane formalism, by 
presenting the action and the field equations. Furthermore, we also study the geodesic equations along the fifth 
dimension, in order to explore the particle motion in the neighbourhood of the brane. In Section [TTT| we present new 
soliton models and discuss their fimdamental properties, for instance, by exploring the broken Z 2 -symmetry character 
of the solutioiLS and the confining effects of the scalar field on the brane. In Section |IV[ we analyse the stability 
of these brane models, where the metric and the scalar field are perturbed about the static brane and the resulting 
equatioiLS of motion are linearized in the proper gauge. Finally, in Section |V| we present our concluding remarks. 


II. THICK BRANE FORMALISM 


We consider a thick brane, embedded in a five-dimensional (5D) bulk spacetime, modelled by the following action 




( 1 ) 
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where is tlie metric and the .scalar cvirvat\ire in the bulk: is a dilatoii held living on the bulk and V[<^) 

is a general potential energy. Note that we are Hiding = SttGs = 2. 

The simplest line element of tlie brane, embedded in the 5D bulk spacetime can be written as [T^ : 

ds\ = 

= dw^ + e^^{dx^-^dy^-\-dz‘- di^), ( 2 ) 


where A is tlie warp factor which depends only on the five-dimensional {5D) coordinate it?. For the scalar held ^ with 
the potential the 5D energy-momentimi tensor is given by 


Tab = dA<^dBy> — 9 ab 


-dc^p d^’p + V{ip) 


where gAB ^ depend only on w. 

The 5D gravitational and scalar held equation.^ take the following forms 


3^"+ 6.1'^ = -Kie-^-^Too =-/c 


+ V{tp) 


,/2 


M' = K^r44 = «:5 
dV{ip) 


1 




+ 4>1V = 




(3) 

(4} 

(5) 

( 6 ) 


respectively, where the prime denotes a derivative with respect to w. 

In order to obtain a hrst-order equation, we introduce an auxiliary fimction W according to (S EZHUJ, which 
demands: 


A' 


wliile V’((^) takes the following fonn [51 1571 - 110] : 


1 dW{p) 

2 dp ' 


VM = I 





(7) 

( 8 } 

(9) 


The Too distribution on the bulk, which will be analysed in detail below, is given by [10] 


Too = e 


2A 


1 

2 



+ V{p) 


( 10 ) 


It can also be shown that, for models with an inhnitely thin brane and Dirac delta distributions, the energy density 
is equal to the cosmological constant of the bulk plus the energy density on the braiie, i.e., e = Af k6{w). 

Moreover, it may be iiLstructive to calculate the geodesic equation along the hfth dimension in a thick brane, in 
order to investigate the particle motion near the brane [H]. As mentioned before, thick brane models considered in 
this paper do not have Dirac delta singularities which enable easier direct calculatioiLS. To this end, we start with the 
geodesic equation: 


wliich leads to 


dx^ 

4 dx'^ dx^ 
dr^ ^ dr dr 


d 


{-2e^’^i) = 0, 


dT 

w + = 0 , 


( 11 ) 


w + = 0 . 


( 12 ) 


where ci is a constant of integration and the function /(u?) is defined as 


f{w) = 


( 13 ) 
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Equation is a second order diiferential equation for w and its sohition depends critically on whether f(^w)/w 
is pasitive or negative. For positive values of f(w)lw one obtains periodic (exponential) solutions, respectively. Note 
that the periodic (negative) motion indicates particle confiiienient near tlie brane, while the exponential solutions 
implies tliat the reference point is imstable. However, tins may point to the possibility that u? = 0 is different from 
the localization of the brane. In a periodic situation, by introducing a new quantity F{w) = c^A' one can 

write the geodesic equation in the following form 

u} + F(u;) = 0. (14) 

The equilibrium point wq satisfies F(ufo)=0. On the other hand, by expanding F(u;) around icq, we have 

F(u;) = F{wo) -f F'(u;o)(tc - u;o) + ..., (15) 

and the geodesic equation lea<ls to 

w + F'(u?Q)(tz; — u^o) = 0. (16) 

Taking into accotint a change of variable w = w — wq, the geodesic equation reduces to 

w + = 0, or w + = 0, (17) 

where Q. = ^/F^wq), provided that ^'(ico) > 0. 

It is essential to emphasise that in the RS-II brane model, the A'A' zero mode corrasponds to the massless graviton 

and the massive modes form a continuum which result in a small correction to Newtonian gravity at large distances 

1321 US]* Free particles are only affected by the gravitational field and not directly by the scalar field. Any field or 

particle which has a direct coupling with the scalar field, will be further affected by extra force from the scalar field. 

2 

Moreover, even an exponential potential like reduces to a harmonic potential for small amplitude oscillations 

as 1 — aw'^ + 0{'w^)). 

In the following section, we explore several models for thick branes and we will employ the linearized geodesic 
equation (flTf for each model. 


III. SOLITON MODELS FOR THE BRANE 
A. Sine-Gordon-based models 


The sine-Gordon (SG) model is a well-known integrable model which ha.s fomid interesting applications in various 
disciplines [S1I35]« Iw fact, single and multiple (topological) soli ton solutions of this system are found analytically 
through different mathematical methods [S1I35]« The self-interaction potential for tills model reads 

^ “ co,s(6ip)] . (18) 

where a and b are free parameters of the model. When considered as the brane potential, liowever, this potential 
should be modified to become consi.stent with the Einstein equations. 

The SG system has the following exact static kink sohition [44| ; 


fp{w) = - arctan , 


(19) 


wliich is plotted in Fig. [^a), for various values of parameters a and b, which correspond to branes with different 
thicknesses. The formalism of onr investigation is to keep the soli ton solution of the fiat space nonlinear equation 
and modify the scalar field potential in such a way that the soliton sohition remaiiLS a sohition of the full gravitating 
system. This is why the soliton solution remaiiLS the same. The form of the potential, however, changes accordingly. 


Taking into account the scalar field given by Eq. (10 1 , and plugging it into the field equations, we obtain the following 
quantities 

\V{<p) = 


exp(2A) = 





. ' U J. 


2 _|_ '' 

02\/a5u; 

le 


,2y/aB\ 


( 20 ) 
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where tlie warp factor is plotted m Fig. [^a). The correspoacliiig modified potential for this model is given by 


TT/ X 2a . . 


64a 

W 


cos 


T 2 


( 21 ) 


wliich is depicted in Fig. [^a). Notice tliat this potential liaii two series of non-degenerate vacnna. as in the DSG 
(double sine-Gordon) system potential |3S]. However, in the limit of & ^ a these vacnua tend to the same value 
(become degenerate), sucli ajs the potentials used in [TUI[5^ . 

This system leads to tlie following Too- from wliicli the energy density can be obtained, 


Too = 


j _j_ \ ^ ^ 

«3v^u; / 




wliich is plotted in Fig. Qa). The energy momentum tensor is calculated accouling to: 


T^u = - g^u 




( 22 ) 


(23) 


where the modified potential is used for V. The metric and the energy momentum tensor are checked to satisfy the 
Einstein equations. Moreover, note that the potential for any soliton model may be shifted by a constant, without 
affecting the soliton solutions. The mininumi value of the potential (i.e., the cla.ssical vacuum), if negative, leads to a 
negative energy density. This can be avoided by a<lding a positive cocustant to the self-interaction potential, when the 
scalar filed is not coupled to gravity (i.e., in fiat spacetime). In curved spacetime, however, tliis constant is non-trivial 
and plays the role of a casmological constant of the bulk. Of course, a positive cosmological constant violates the 
strong energy condition. 

Note that the energy density is localized at the brane and the thickness of the latter is given by 


A = 



In this model, tlie Ricci and Kretschmann scalars are given by 



and 


A' = 


16384 [>-60 

“”81 



(24) 


(25) 


(26) 


respectively. It can be seen that there is no singularity in the Ricci scalar and/or Kretschmann scalar. In the limits 
of ID —> ± 00 , the Ricci scalar becomes 


lim i? 

u;—►-hco 

lim i? 

ttf—V —CO 


0 , 

5120 a 


and tlie limit of u? ^ 0 yields 


lim R = 


64(-20 + 362)a 
9 63 


Moreover, the mixed Einstein tensor components are given by: 


-^0 _ /^l _ _ /^3 _ 

Uq — — U 2 — U 3 — 


32 a (-16+ 3^2'^'62^ 



(27) 

(28) 


(29) 


( 30 ) 
















a 


6 


Gi = 


512 




J q2V 


2 ’ 


respectively. Note that all the components of the EiiLSteiii tensor in the limits ti; —>• ±oo become 


liin Gi 

w—*+oo 

liin Ci 

W—i — OO 


512 a 


(31) 


(32) 

(33) 


However, in the limit ofn;—►O (for;i = i/ = 0,l,2,3) the EiiLStein tensor is given by 








(34) 


and one can interpret it as the cosmological constant on the brane, i.e., Gy oc A(5j, with A = 

These results can be interpreted in that we have a broken Z 2 -syimuetry in the bulk, ?ls the two sides of the brane 
differ completely. On the right (li) ^ the Einstein tensor and consequently the cosmological constant of the 

bulk vanish, so the bulk is asymptotically Minkowski. However, on the other side of the brane, these quantities are 
nonzero and equal to the constant value 512a/(36^), and as a result the bulk would be de Sitter. The Ricci scalar and 
the Einstein tensor component Gq are plotted Fig. [^a) and Fig. [(^a), respectively. 

Furthermore, by calculating the field equations, 

Gab = kIT'ab- (35) 

one verities that = 2, which is consistent with the iLsual normalization notation |^|^|2B]. As pointed out in 
DeWolfe et ai 0, in the stiff limit where ab cc, the wall reduces to the step function and the energy density 
approaches a J-function (see Figs, [^a) andj^a)). 

The local confining gravitational field of the brane is best observed by looking at the geodesic equation of a test 
particle moving only in the direction of the extra dimeiLsion, a.s explained in the previous section: 


.8 32^2^0(3^2 igs .8 32^23i^v^ 


(36) 


Tins proves the confining effect of the scalar field, and ci is an integration constant [see Eq. 12 . For small amplitude 
oscillations, the relativistic motion reduces to a classical motion in a Newtonian classical potential. Since we have 
considered the potential up to second order, relativistic effects can be ignored and the corresponding quantum energy 
levels are therefore those of a non-relativistic quantum particle. So, if interpreted as a quantum oscillator, one can 
assign an energy to each quantum state given by En = (n + l/2)hiAj. where 


I. n rE7^—^ ,/832^2si^a(362 + 16) 

Huj = n = y/F'iwo} = tiy --p- 

However, an important point is in order. In what follows, we explore the confining effect of the brane, only up to 
second order in the potential. Even if the classical test particle is confined up to this order, large amplitude oscillations 
will involve nonlinear effects and this might exploit confinement. Quantmn mechanically, the full nonlinear potential 
might lead to timneling and thus de-confinement. It is well known that massive KK modes are not confined to the 
brane. 


(37) 


B. cp"*-based model 

The models are well known for having simple soliton-like solutions, although it is not strictly integrable like the 
SG system. This model is also the central ingredient in the Goldtone and Higgs mechanisms. Spontaneous breakdown 
of the Z 2 symmetry in the complex version of the model leads to the appearance of the Goldstone mode and once 
coupled with a Gauge field, it causes the Gauge boson to acquire mass 02] • This model is therefore frequently used 
for building thick branes. 

For the -based model, we have JE] 




( 38 ) 
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where a aiicl ^ axe coiLSttuits. The kink sohitioii reads 

^{w) = athnh{0w), (39) 

wMch is depicted in Fig. [^b). Proceeding in a similar inamier as in the previous case, we have the following solutions: 


O Ct 

A = -io2i,i[cosh(^iii)] + i«2 — J 

9 9 cosh (pn;) 


exp (2.4) = [cosh(^ii?)] exp 


2^2 1 
.9 cosli^(/iii;). 


respectively, and the potential is obtained as 

1 


( 1 - r_ ) _ { 3 - 


9 \ 2 . 

‘f \ ^ .2 2,02 In 


a* 


a 


27 


(40) 

(41) 

(42) 

(43) 


wliich is plotted in Fig. |j^b). It is seen that while V{<^) was 0{<p^), V(^) is 0(^®). 
The corresponding Tqq is given by 


1 exp 1 — 5 a^tanh‘‘( 0 ti; 


2^.2 


^00 = — 


27 




— 4rt^ — 27 cosh^(£^ii;) 


coslr(/3n;) — 12 fV^ coslr(/3n;] 


(44) 


wliich is depicted in Fig. ^b). where one verifies that the energy density is localised at the brane. The brane thickness 
becomes A = 

Moreover, one can show that the Ricci and Kretschinann scalars are given by 


R=- 


16 a 


2o2 


81 cosh® 

wliich is depicted in Fig. [^b), and 


— ISrt^ cosh^(5uf) — 5 q^ + 20fv^ cosh®(/?nr) — 27 cosh^(^tij)j, 


(45) 


4 rj4 


K = 


64 


QOa"^ co^\i^{lSw) — cosh^(^ii?) — 240('t^ cosh*(/?n; 


6561 cosh^^(3ii?) 

+160^** cosh^^(^nj) + 60(v'^ cosh^(^n;) + lOo^ + 324a^ cosh^(5nj] 
+108rt^ cosh^(0tn) - 432o^ cosh®{^nj) + 729cosh‘^(£iii;)j, 

respectively. Tlie mixed Einstein tensor components take the following form 

2 


(46) 


G o r^l _ /-t3 

0 — Ltj — '-r2 “ CJg — 


[ — 4a^ - 12o^ coslr{^nj) 


^ 27cosh®{^nj) 

+16a^ cosh®(^ii?) — 27 cosh^(3ti;)], 


(47) 


Gi = 


8 [—3cosh^(5m) — 1 H- 4cosh®(3uf)] 


27 


cosir (/^u;' 


(48) 


where the Go component is depicted in Fig. [^b). 

As for the previous SG model, we determine the limits it? —V ±ao for all the components of the EiiLstein tensor 
components, which are given by: 


li- 

u)^+oo 2i 

lini Gi = 

ui—♦—oo 2i 


(49) 


( 50 ) 
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FIG. 1: Solitori solutions as a function of the five-dimensional coordinate w for the following models; (a) SG for 6 = 1, (b) yj* 
for a = 1 and (c) y?® for a = 1 systems. Dashed, dotted-dashed, and continuous curves correspond to solitons with decreasing 
branc thickness. In the limit of an infinite a/Q parameter, the soliton approaches the step function. 
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FIG. 2; The plots depict the warp factor as a function of the fifth dimension for the (a) SG with a = 100 and 6 = 1 , (b) 
with a = I and ^ = 100 and (c) with a = 1 and 5 = 100 systems, respectively. For the SG and systems there is an 
asymmetry between the two sides of the brane and the Z 2 symmetry is broken. 


and ill the limit of ^ 0 the Einstein teiLsor components (for = i/ = 0,1,2,3) takes the form 


lim Gj; = - 20 ^ 0 ^ 6 }^. 

w —*0 


(51) 


Note that in this model the cosmological constant on the brane A would be Taking into account all of the 

above considerations, we verify that the Einstein equatioiLS are given consistently by Gab = where K 3 = 2. 

The geodesic equation for a test particle moving in the direction of the fifth dimension one obtains 




(52) 


wliich corresponds to a linearized quantum mode of energy hu) = il = 


y/F'{wo) = ■^CiaQ. 


It is seen tliat the resulting potential for tlie model is an odd fmiction witli respect to u;, while the A fmiction 
is even. This property is not verified in the SG and models, where the latter is discussed below. This results in 
an nsyininetry in the corresponding properties (such as the energy density). However, by an appropriate selection of 
the model parameters, one can restore the Z 2 symmetry in the and cases, which is commonly used in brane 
models. 
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FIG. 3: The plots depict the modified soliton potential as a function of the scalar field for (a) SG with a = 6 = 10, (b) 
with Q = 0 = 1 and (c) with a = 0 = I systems. The potential of the SG aud systems have non-degenerate vacuua. In 
contrast, the ip* potential has degenerate vacuua and this leads to a stable, topological solitonic brane. For and p^ systems 
the number of extrema in the Figure is only a result of the range of the plot. In a wider range plot, other extrema appear: 
these may render the system globally imstable. We have only claimed local stability. 



FIG. 4: as a function of the fifth dimension coordinate u; for the following systems: (a) SG with b = 1, (b) with a = 1 
and (c) with a = 1. The dashed, dotted-dashed and continuo\is curves correspond to increasingly thin branes. The brane 
becomes infinitely thin (delta function) as the parameter (i/0 tends to infinity. 
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FIG. 5: Ricci scalar as a function of the fifth dimension coordinate w for the following systems; (a) SG with 6 = 1, (b) p* with 
a = I and (c) p^ with a = 1. The Ricci scalar is different on the two sides of the brane, due to the effect of the warp factor. 
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FIG. 6: The plots depict the Einstein tensor compouent Go for the (a) SG with a = b = 1, (b) with q = 5 = 1 and (c) 
with ct = 0 = 1 systems. Note that this quantity approaches different constant values for the SG and systems, while the 
system is Za‘Symmetric. 


C. (^^-bascd model 


For tills model, we have the following potential: 

= (53} 

and aiS a result, the kink solution is given by 3H] 

ip{w) = ^ , (54) 

^Qi-y/^Q0w) 

where a and 0 are constant (as in the model). The kink solution is depicted in Fig. [^c). For this model the 
potential W, and the warp factor A are given by the following expressions 


W(^) = 


\/2 — (p‘ 


a 


A = - 


12 


1 J_ 


In 


1 4- e 




,—*/2ai0w 


exp(2^} = 


,—Wy/2da 




g 6 


respectively. Thus, for the system the self-interaction potential takes the form 

,,, , 1 

= 4-^- 

wliich is depicted in Fig. [^c). Note that V((^) lias raised to 0(^5®). 

The energy-momentum teiLsor coniponeut Too is given by: 


^00 = —— 


1 / Q—w^/2Qa \ ^ g ® i+«-wVTflt 




24 1 + j g-u.v/2S^^ 

+ 4a2e-“’'^^“ + j, 

wMch is plotted in Fig. Qc). The thickness of this hrane is given by A = {\/2a^)~^. 


(55) 

(56) 

(57) 


(58) 


( 59 ) 
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For tlii^ system the Ricci scalar is given by 


1 ^4/32«ttfN/2Sa r _ 

_ i- a P e \r,n^2^wV2Bc 


R = -^ ~ ^ ' -^|20Q'e“'''-=°“ 


+ l) 

+50^ + 20a2 _ 48 ^ 

wliich is plotted in Figs, [^e), and the Kretsclmiann scalar by 

j ^ ^2wy/2^< 




1152 + 5o^ ('e6u.v/2Sc.'\ 4_ 4 q^4 / 

('e..v/2Sc + i)’L V 7 V ; 

+ 120rt-‘ + 160o^ (e3..v/2^)a^ 8 q^4 ^^2wV20^-j 

-480rt^ 1152 ^e2u,v/2d=i^ _ 7g8^2 ^g2u.v/2Sa^ 

-384«2e“'^^“ - 96f»^ ^e4,^^/2^a^ + 2304e“'’^^“j. 
The mixed Einstein tensor components are given by 


C o /~<1 /^2 /~<3 

0 — Ui — U2 — “ 




{a 


2„3wV2^a 


12 + l) 

- I2e“’'^^“ - I 2 ), 


(60) 


(61) 


(62) 


Gt = 


^ a^fi‘2^'2wy/20a ^2 wV20q _|_ 


12 


(^QWy/20a ^ 


(63) 


where the Cq component is depicted in Fig. [^c). These tensor components reduce to the following in the limit of 
w —> ± 00 : 


Ihn Gi = 

U)^+oo lZ 

lim Gg = 0, 

tL—♦—OC 


and in the limit o? w ^ 0, the Einstein tensor is (for /t = i/ = 0.1.2,3): 




As in the previo\is cases, the Einstein equation in the bulk Gab = ^I^ab ^ found self-consLstently with = 2, 


(64) 


(65) 


and which leads to the following geodesic equation 

.. . 1 


tc + ct — exp 
*192 ^ 


12 


(21n2 + l) 


1 


-^exp 


(21n2 + l) 


(S + 3«2) a) 
a^0\/2. 


( 66 ) 


The quantum mode energy is thus given by 

fk*j = i} = y/F'{wQ) 


Cr^VzJexp Cg (2lii(2) + 1)) (8 + 3a2), 


( 67 ) 
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IV. BRANE STABILITY 


III this .section, we examine small perturbations about the soliton branes obtained in the previous sections. To this 
end, the metric and the scalar field are perturbed about the static braiie and the resulting equations of motion are 
linearized in the proper gauge The linearized equation tinns out to be a Sclirddinger-like equation with 

a potential U{z) which determines the linear modes. Unfortunately, for the models considered, this potential is too 
complicated to be handled analytically, or to be used for finchng the corresponding modes. Therefore, we will only 
examine the potential near its mini mum up to second order in 2 . 

In order to check the stability of the brane models, we iLse the standard scalar-tensor-vector (STV) decompositions 
p2] . For this purpose, a new variable 2 is coiLsidered in such a way that dz = so the metric can be 

considered as 1521: 


9AfN = 


( 68 ) 


where a(z) = On the other hand, one can linearize the action of the system given by Eq.Q along with the 

metric, via and iJpA/iv = which is given by [52]: 


^ I d^xa^ 


dMliNpa^’h’'’^ -d^'hd^h 


A/JV 


-IrdphMNd^h^’'^ + l-d'’hdph + 3- 

+ 2 K 5 j. 


(69) 


where On, 9^ = the prime denotes a derivative with respect to 2 . Moreover, 

by introducing the following vector and tensor perturbatioiLS 


70) 


— d^F-\-G^, 

d^duB , 

where (7^ and are the transverse vector and tensor perturbations, respectively, one obtaiiLS 

= 0 = 

=0=1?;;. 

By using this STV decomposition method, one can decompose into the following decoupled parts 


:7i) 


<5'2ls, 






^ a3 


where 

The scalar perturbations of the action, in turn, lead to two parts |521155) : 


(72) 


73) 


:74) 


with ijj = F — and 




jca/ar—! 


r-2 ~ 2 J ~ + 6l^V* 


-3^/i,,(/i', + 4^') + 24 + a^C^^[5^)' 

+2/t„v’'V + + 4vp')-5v - (V)') 


(75) 
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After eliminating hi action by taking into accoiuit Eq.<j74|f. and doing some simplifications, one obtains 

HUSH]: 


= I rf'':re|n'‘’'e + d"-ye|, 


(76) 


where G and 0 are gauge invariant variables and specific fmictions of o, and are given by 

-2 




.. 3 a 

0 = 


respectively. Finally, the equation of normal modes of the linear perturbations are: 

vector: 


(77) 


teiLSor: 

scalar: 


+ d;, - = 0. 

a'2 

a(*)C + S''-^S = Q. 

u 


(78) 


Since the vector perturbations only have zero modes, solutions are stable against vector perturbations 32 • However, 
for tensor and scalar modes another decomposition should be introduced 321) namely, 


^{x\z) = e'9'"'€>,(2). 

where is the TT polarization tensor and Pp(2) and ^q{z) satisfy the following equations: 

AtA\pp = mlpp. 


(79) 

(80) 


with nip = —p^p^, and 


A,Al^g = 


X = A + 

“ T ^ a ’ 

dz as 
, d 0" 

= :r + -F* 
dz 9 


Equation (j^ is a Sclirddinger-like equation, and takes the form 

- P‘p + W(2)pp = mlpp. 


(81) 


(82) 


(83) 


where the effective potential VV( 2 ) is given 


W(2) 







(84) 


The effective fo\ir-dimensional gravity is determined by the spectrum of the tensor A'A' inodes pp 32* instance, 
for zero mode po with oio = 0 and a normalizable po leads to the fom-dimensional Newton’s law [T^150] . 

On the other hand, in order to study the stability of the branes, some authors choose an “axial gauge” where the 
metric is perturbed as p iSTIlin] : 


da^ = eh^i,)dx^dx^ — dw^. 


(S5) 
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FIG. 7: The coefficient of the term in the potential of the linearized Sclirddinger equation as a function of the free parameters 
Q and b for the SG system. Negative values correspond to a first order instability. 



FIG. 8: The coefficient of the term in the potential of the linearized Sclirddinger equation as a function of the free parameters 
a and for the system. The coefficient is everywhere positive, signaling linear stability. 


Here represents the background metric, denotes the metric perturbations, and e is a small parameter [35]. 
They also consider the transformation where = (p(xyw) ^Oj. Moreover, in order to render the metric 

conformally flat, one can choose dz = In this case, the corresponding Sclirddinger equation takes the form 
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FIG. 9: The coefficient of the a" term m the potential of the linearized Sclirddinger equation as a function of the free parameters 
a and 8 for the system, Negative values correspond to a first order instability. There are also vast patches in the parameter 
space which have almost neutral stability. 


laiisiEzniniEii: 

- ( 86 ) 

where the potential is given by: 

U{z)^-\k+\^^ + \a". (87) 

As one can check that, this potential and the V’ function are the same af? W( 2 ) and p in Eq. Note that A is 

a cosmological constant on the brane, which could be positive, negative or zero corresponding to the 4D spacetime 
being do Sitter (^ 64 ), anti-de Sitter (AdS 4 ) or Minkowski {A/ 4 ) [33 [35|. Besides, it is notable that the Hamiltonian 
corresponding to Eq. can be written in the form (251 IHZl !SI]* 


wliich is obvioiLsIy Hermitian and therefore leads to real k {k^ > 0). Accordingly, there are no unstable tachyonic 
excitations in the system [iUl 151) . 

The solution for the zero modes (fc = 0) is pnil421[5T| ; 

\p[z) = Ne~^ (89) 


where iV is a normalization factor [iHl 1121 [ST] and satisfies: 


1 


/ 


+00 


rf-z|»/'o(z)|' 





jV2 r+°° 

~ 7-00 




(90) 


where y = Iw is & dimensionless variable. The a.symptotic behavior of the solutioiLS at large y are checked and the 
result is that only the system ha.s a normalizable zero mode and thus stable. On the other hand, aa in quantum 
mechanical systems, we may check for the stability of the system via the existence of a real frequency, in bound states. 
Since the potentials for the three systems considered in this paper, are too complicated to be solved analytically, we 
found the corresponding ground state eigenvalues via expansions in terms of the fifth coordinate w. One can dedtice 
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the stabiUty up to 0(w^) by looking at the .sign of the iv'^ term. It is seen in Figures tJo that this coefficient is 


everywhere pasitive (stability) for the system, while there are regions of the parame ter space where the coefficient 
is negative (instability) for the SC anil systems. However, as noted in section III Al this conclusion is not decisive, 
since higher order effects might have drastic effects. The difference between the results of these two approaches is 
probably caused by the inevitable approxiniation.s iLsed in the analysis. 


V. CONCLUSION 

In this work, we obtained exact thick brane models iiLSpired by well-known nonlinear systems, namely, the sine- 
G or don (SG), and models. The confining effect of the scalar field in all these three models were confirmed by 
examining the geodesic equation for a test particle moving normal to the brane. In particular, it turns out that the 
modified potential for the SG system resembles tliat of the double sine-Gordon {DSC) system, wliile those of and 
became and respectively. We have extended previous brane models [El[131122] ba.sed on SG and solitons 
taking into accoimt different parametrizations. The similarity of the model with the generic Higgs model makes 
this choice particularly interesting, especially as the resulting potential is an odd function of the fifth coordinate and 
the Z 2 symmetry is respected. We liave studied the model for the first time. TliLs model is interesting by its own 
right, since unlike the model, we have two pairs of solitons and anti-solitons which live in different topological 
sectors. 

In the cajse of the SG model, we have used a more general form of the potential compared to the one used in (10] 
and [33|. The resulting brane does not have Z 2 symmetry, in general, where the center of the brane may be displaced 
from n; = 0 and the potential will not be an odd function of w in general. However, by a suitable choice of the model 
parameters it is possible to make the vacuua of the effective potential degenerate, in which ca^se the Z 2 .symmetry is 
restored. In the case of the model, however, we could not restore this syimnetry via re-parametrkatioii. Finally, 
using standard procediues, we examined the stabihty of tlie thick branes, by determining tlie sign of the tenn in 
the expansion of the potential for the resulting Schrodinger-like equation. It tiuiLS out that the hrane is .stable, 
wliile there are unstable modes for certain ranges of the model parameters in tlie SG and branes. 

We considered tlie limiting case in which the brane tends to zero thickness and approaches a tliiii brane. It should 
be noted that the topological stability of the soliton brane remains valid even in this limit (at least at tlie classical 
level). An interesting question would be whether the thick brane continuoiLs metric develop.s a discontinuity and 
whether the Israel junction conihtioiis will be satisfied in this hmit. Although one would expect intuitively that this 
is tlie case, we have not worked out tlie detailed calculations. This issue will be explored in separate paper. 
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